Effect of Defects on the Line shape of Electron Paramagnetic 
Resonance Signals from the Single-Molecule Magnet Mni2: A 

Theoretical Study 

Kyungwha Park^'^'*, N. S. Dalal^'t, and P. A. Rikvold^'^'^ 
^School of Computational Science and Information Technology, Florida State University, 

Tallahassee, Florida 32306-4120 
'^Department of Chemistry, Florida State University, Tallahassee, Florida 32306-4390 
^ Center for Materials Research and Technology and Department of Physics, Florida State 

University, Tallahassee, Florida 32306-4350 
(February 1, 2008) 



1 



Abstract 



We herein estimate the effect of lattice defects on the hne shape of electron 
paramagnetic resonance (EPR) signals from a single crystal of the S = 10 
single-molecule magnet Mni2 with the external magnetic field along the crys- 
tal c axis. A second-order perturbation treatment of an effective single-spin 
Hamiltonian indicates that a small, random, static misorientation of the mag- 
netic symmetry axes in a crystalline lattice can lead to asymmetric EPR 
peaks. Full spectra are simulated by calculating probability-distribution func- 
tions for the resonant fields, employing distributions in the tilt angle of the 
easy axis from the c axis, in the uniaxial anisotropy parameter, and in the g- 
factor. We discuss conditions under which the asymmetry in the EPR spectra 
becomes prominent. The direction and magnitude of the asymmetry provide 
information on the specific energy levels involved with the EPR transition, 

the EPR frequency, and the distribution in the tilt angle. 
PACS numbers: 75.50.Xx, 76.30.-v, 75.45.+j, 61.72.Hh 
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I. INTRODUCTION 



The nanomagnetic compound [Miii20i2(CH3COO) i6(H20)4]-2CH3COOH4H20 (abbre- 
viated hereafter as Mni2)/ first synthesized by Lis/ consists of a core of twelve Mn ions, 
hnked together by oxygen atoms from the bridging acetate hgands. The Mn ions form the 
magnetic core, and are shown schematically in Fig. 1. The Mn ions are of mixed valence: 
eight of them (Mn^+) with four unpaired electrons {S — 2) form the outer ring, and are 
ferromagnetically coupled within themselves, to lead to a total spin of S*! = 8 x 2 = 16. The 
other four ions (Mn^+) each with S — 3/2 form the center of the core and are arranged in a 
Mn-0 cubane type structure. These are also ferromagnetically coupled within themselves, 
with a total spin of 5*2 = 4 x (3/2) = 6. The total spin on a Mni2 molecule has been deter- 
mined to be 10, which can be understood as arising from a strong antiferromagnetic coupling 
between the two Si and 5*2 sub-systems. The whole molecule may thus be considered as 
a 5* = 10 particle as far as the ground state is concerned, which is of primary interest in 
the present work. The acetate ligands are covalently bonded via the oxygens to the Mn 
ions. There are two acetic acid molecules and four H2 molecules solvated in the lattice. The 
magnetic moments of the Mni2 molecules are most easily aligned along the crystal c axis, 
which is thus the easy axis for a perfect crystal (without imperfections/dislocations). There 
is a zero-field energy barrier of approximately 65 K against magnetization reversal. ^'^ The 
effective distance between different Mni2 molecules is approximately 14 A.^'^^'^'' 

Earlier, we have examined the effect of possible defects on the line widths of electron 
paramagnetic resonance (EPR) spectra for the single-molecule magnet Mni2,^^'^^ and have 
compared our calculated results with experiments. The imperfections in large complex 
molecules such as Mni2 could be of many types, such as dislocations^ and the disorder 
caused by the various orientations of the acetate ligand and of the solvent molecules in 
a unit cell.* So far the origin of the possible defects has not yet been clarified. In the 
study, ^^'^^ we modeled that the defects could be represented by the distributions in the 
uniaxial anisotropy parameter D (second-order anisotropy) and the gi-factor. To calculate 
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the EPR linewidths, the intermolecular dipolar interactions and the distributions in D and 
g were considered. For Mni2 the intermolecular dipolar interactions (~ a couple of hundred 
gauss) are substantially smaller than the the effect of the distributions in D and gi (~ a 
thousand gauss). The detailed estimates of the dipolar interactions and the effect of the 
distributions in D and were given in Refs.^^'^^ as a function of temperature, resonance 
frequency, and the energy levels relevant to the EPR transition. Our calculated results^^'^^ 
agree well with experiments. ^^'^^ 

In the present study, we describe our theoretical investigation of the effects of lattice 
defects on the signal line shape of EPR spectra for Mni2. We propose that the lattice 
defects may cause a small, static, random misorientation of the magnetic symmetry axes 
of each molecule in the lattice. Here we examine the full microwave absorption line shape 
as a function of resonance frequency and the energy levels relevant to the EPR transitions, 
including Gaussian distributions in the tilt-angle (the angle between the magnetization easy 
axis and the crystal c axis) as well as the distributions in D and g. Then we find that 
consideration of the Gaussian distribution in the tilt angle leads to a line-shape asymmetry 
in the spectra. Conversely, it might be possible to learn about the defect structure from this 
asymmetry. 

This study was prompted by three factors. First, Mni2 and its analogs, with their 
high-spin ground states of 6" = 10, exhibit novel properties, such as macroscopic quantum 
tunneling (MQT),^'^^ whose origin has not yet been fully understood. A new model of MQT, 
especially for Mni2, involves a direct role for lattice defects.^ Second, Mni2 has recently been 
proposed as a material for potential use in quantum computation.^^ For this application, a 
detailed knowledge of the hne shape of the EPR absorption spectrum would be required. 
The proposed method can be applied to microwave transitions between the various spin 
levels. Third, many earlier high-field EPR studies of Mni2 and related systems have reported 
asymmetric line shapes, but the origin of the asymmetry has not been clarified. ^'^^ It was 
assumed that the asymmetry might be an instrumental artifact, resulting from the existence 
of standing waves in the waveguides. ^'^^ 
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Asymmetric EPR line shapes have been known since the early days of EPR spectroscopy. 
For example, Bleaney and Rubins considered forbidden hyperfine transitions and noted that 
EPR line shapes from powder samples, which have a distribution of particle orientations, 
can be asymmetric.^^ Although we use a standard spin Hamiltonian [Eq. (1)] and routine 
perturbation technique, our study is distinct from such earlier studies in the following ways: 
(a) We consider single-crystal samples, with measurements along a well-defined direction 
such as the crystal c axis, (b) Our system involves the existence of strong intramolecular spin- 
exchange effects within the Mni2 molecule (which is an order of 100 K),^ which obliterates 
the hyperfine structure and the effects of any forbidden transitions, (c) The system and 
experimental situation are characterized by a dominant uniaxial anisotropy term in the spin 
Hamiltonian and a strong external magnetic field, and are thus quite applicable to high-field 
EPR. (d) We examine systematically the asymmetry in the EPR line shape as a function of 
the energy levels between which the microwave transition occurs, the microwave frequency, 
and the distribution in the tilt angle. Future experimental tests of our predictions for the 
EPR line shapes of single-molecule magnets should be fruitful. 

The remainder of this paper is organized as follows. In Sec. II the model is described, and 
in Sec. Ill the resonant field is calculated as a function of the rotation angle 9 (which is the 
tilt angle of the molecular easy axis with respect to the crystal c axis), and the direction and 
the magnitude of the asymmetry in the resonant fields are discussed. In Sec. IV we calculate 
the probability-distribution functions of the resonant fields, including the distributions in 9, 
D, and g, to predict the asymmetry in the spectra. In Sec. V we present a discussion and 
our conclusions. 

II. MODEL 

In the presence of defects or impurities in a sample, Mni2 molecules can be displaced or 
rotated from their normal positions. Consequently, each molecule sees a slightly different 
crystal field caused by the surrounding molecules, compared to that seen in a perfect crystal. 
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We propose that this shghtly different crystal field seen by a molecule results in a small, 
static, random misorientation (rotation) of the magnetic symmetry axes of the molecule in 
the crystalline lattice. The degree of the misorientation can be quantified as the magneti- 
zation easy axis of the molecule is rotated by an angle 6 away from the crystal c axis. The 
majority of the molecules are assumed to have their easy axes almost aligned with the c 
axis, so the tilt angle 9 is assumed to have a Gaussian distribution about zero with a small 
standard deviation. Hereafter a, 6, and c denote the crystal axes, while x, y, and z denote 
the molecular (magnetic) anisotropy axes of a single molecule. In this study, we consider 
the case of varying the magnetic field at a fixed EPR frequency (roughly 50 to 500 GHz) 
when the external magnetic field is applied along the c axis. 

We confine our study to the uniaxial molecular magnet, Mni2, because of its very 
small transverse anisotropy compared to its uniaxial anisotropy. Additional transverse 
anisotropy terms or small transverse magnetic fields at ^ = may also produce asymmetry 
in the spectra, because of nonlinear relationships between the energy cost of an EPR tran- 
sition and the sweeping field. For Mni2, the negligibly small transverse anisotropy enables 
us to concentrate on the asymmetry effect caused by a distribution in 9 only. 

Our goal is to investigate how static, random rotations of the magnetic anisotropy axes 
caused by defects affect the line shapes of EPR spectra as functions of frequency and energy 
level, and to predict EPR spectra which may be compared with experimental data in near 
future. For the single crystal of Mni2 examined in Refs. 13, 18, it was found that the effects 
of D-strain and (yf-strain are more significant than dipolar interactions. ^^'^^ Therefore we need 
to combine the effects of I?-strain and g^-strain with the effect of static, random rotations of 
the easy axes in order to obtain realistic spectra. After analytically calculating the resonant 
field as a function of 9, D, and g, we numerically obtain the probability distribution function 
of the resonant field, assuming Gaussian distributions in 9, D and g. We do not take into 
account the effects of natural linewidths, of dipolar interactions between molecules, or of 
temperature on the line shapes. 

When an external magnetic field is applied along the c axis, and the magnetic anisotropy 
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easy axis of a single molecule (the z axis) is tilted by Q away from the c axis, the single- 
spin ground-state Hamiltonian is, to lowest order, in terms of the spin operators along the 
molecular axes. 



where D — 0.55 fc^,^ the isotropic g — 1.94,^ \xb is the Bohr magneton, and B is the 
magnitude of the external field. Here, is the third of the standard Euler angles (6*, 0, as 
defined in Ref. 19. (For the rotation matrix, see Appendix A.) There is no constraint on the 
magnitude of V' because it does not affect the energy eigenvalues of the spin Hamiltonian. 

For simplicity, we ignore small fourth-order anisotropy terms^'^ in the single-spin Hamil- 
tonian. This is justifiable because: (i) these terms make rather small energy contributions to 
the eigenvalues, and (ii) some experimental data indicate that the magnitude of the trans- 
verse fourth-order terms for Mni2 may be too small to be measured via, for example, neutron 
scattering.^^'^^ 

We can also write the single-spin Hamiltonian in terms of spin operators along the crystal 
axes. Both forms of the spin Hamiltonian represent the same system. (See Appendix A.) 



To calculate perturbatively, to lowest order, the resonant fields with the spin Hamiltonian 
(1), we assume that the tilt angle 6 is very small. This assumption is probably valid because 
a large misalignment between the crystal c axis and the magnetic easy axis of a molecule 
(the z axis) is not expected for Mni2, since the reported X-ray data^^ on bond angles 
and distances are quite precise, with no mention of mosaicity issues. Then taking V = 
—q/ibB sm9{smip Sx + cosip Sy) as a small perturbation to the rest of the terms in H, we 
obtain the second-order perturbed energy of the level Ms, where Mg is an eigenvalue of 
the spin operator S/. 



n 



DSl — gfisB cos9 Sz — g/isB sin 9 {simjj Sx + cos^ Sy) , 



(1) 



III. RESONANT FIELDS 



S 



(2) _ v\S-AQ{S + Ms + l) ^ ri\S + Ms){S-Ms + l) 



(2) 
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where rj = (giiBB sin 6) /2, and — —DM^ — gjisB cos 9 Mg is the unperturbed energy 
eigenvalue of the level Mg. With the selection rule that EPR transitions occur between 
adjacent energy levels only, we have at resonance hv — ^m]±i + ^m!±i ~ ^m! ~ ^m!' where v 
is the frequency of the applied microwave radiation in the EPR experiment. The resonance 
equation is quartic in so its solutions are too complicated to provide useful information. 
To obtain a simple analytic expression for the resonant field, we substitute the unperturbed 
expression for the resonant field Bq [for example, for the transition between the levels M, 
and M,±l, Bq = {^hu - D(2M, ± l))/^/iBCOs^], for B in Eq. (2). Then we obtain the 
approximate resonant fields as follows: 

BUMs - - l)=So(M, - - 1) - 9^BBisin^eX{.,M,) ^ ^ 

t: cos U 

BUMs - + l)=5o(M. + 1) + gi^BBl^^^Xi.^-Mg) ^ ^ 

v( 'KS-Ms + l){S + Ms) {S + M,-l){S-M, + 2) 
A [y, Ms) — 



hv hu - 2D 

(,5-M,)(5 + M, + l) 



(5) 



hv^2D 

Notice that the sign of the term proportional to X for the transition Mg Mg — 1 is 
opposite to that for Mg — > M, + 1. The first type of transitions, M, — > — 1, represent 
transitions among energy levels in the lower-energy potential well (see Fig. 2). The second 
type, Mg —>■ Mg + 1, corresponds to transitions among energy levels in the higher-energy 
potential well (see Fig. 2). At frequencies approximately lower than 100 GHz both types 
of transitions are observed, while at higher frequencies only transitions in the lower energy 
are observed. This is because at high frequencies the EPR excitation energy exceeds the 
difference in adjacent energy levels relevant to the transitions. 

Our approximation [Eqs. (3) and (4)] is valid when the absolute value of the ratio be- 
tween the second-order and the zero-order terms becomes very small compared to unity: 
\Y{u,Mg,e)\ = IgneBosm'^e X(i/, :fMJ/(4cos^)| < 1 for Mg ^ Mg ± 1. Figure 3 shows 
y(i/, Mg, 9) as a function of frequency v and at ^ = 0.1 rad 5.7°) for the two types 
of transitions. For the transitions Mg — > Mg — 1 (we consider the case that the resonant 



field is positive at ^ = 0, that is, hu — D{2Ms — 1) > 0), the approximation becomes good 
for all Mg in the range of frequencies from 50 GHz to 500 GHz, except for higher energy 
levels than Ms = 1 [Fig. 3(a)]. The approximation is better for higher frequencies and 
larger Mg (lower energy levels). For the transitions Mg — > Mg + 1 (we consider the case 
that D(2Mg — 1) — hu > 0), the approximation becomes better at higher frequencies or for 
5ma//er magnitude of Mg (higher energy levels) [see Fig. 3(b)]. As shown in Fig. 4, the 
resonant fields from this approximation agree well with the exact diagonalization results at 
u = 66.135 GHz.24 

Next we examine the consequences of the nonzero 9 on the resonant fields, Eqs. (3) and 
(4), as functions of the frequency and the energy level. Gaussian distributions in D and 
g provide symmetric distributions in the resonant fields. However, a Gaussian distribution 
in 9 lets the molecules have smaller or larger resonant fields than that for ^ = 0, in an 
asymmetric fashion. This yields asymmetric tails in the average line shapes of EPR spectra. 

To investigate the asymmetry effect as a function of v and Mg, we define the quantity 
A{i^, Mg, 9) = B^^{9 7^ 0) - B^^{9 = 0). As shown in Figs. 5(a) and 6(a), for Mg^ Mg- 1, 
the sign of A{v, Mg, 9 — 5.7°) is positive so the resonant fields become minima at ^ = 
for the whole examined range of frequencies and for all energy levels in the lower-energy 
potential well. Therefore, the spectra will have asymmetric tails in the direction of increasing 
field. The magnitude ol A{y,Mg, 9 — 5.7°) increases with decreasing Mg at low frequencies 
(below about 200 GHz) [Fig. 5(a)] because of high resonant fields, and decreases weakly 
with decreasing Mg at high frequencies [Fig. 6(a)]. For a particular Mg, the magnitude of 
A{u, Mg, 9 — 5.7°) increases with increasing frequency, except for = 3 — > 2, = 2 — > 1, 
and M5 = 1 — > [Fig. 6(a)]. For = 3 — > 2 and Mg = 2 — > 1, the magnitude increases 
with frequency for low frequencies, and then starts to decrease at about 100 GHz and 250 
GHz, respectively. For Mg = 1^0 the magnitude decreases monotonically with increasing 
frequency. As shown in Figs. 5(b) and 6(b), for Mg^ Mg + 1, the sign of A{i', Mg,9 = 5.7°) 
is negative so the resonant fields attain their maxima at ^ = for the whole examined 
range of frequencies and Mg in the higher-energy potential well. Thus, the spectra show 
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asymmetric tails in the direction of decreasing field. The magnitude of M^, 9 — 5.7°) 
increases with increasing magnitude of Mg [Fig. 5 (b)] and with decreasing frequency [Fig. 
6(b)], because the resonant fields increase. It should be recalled that our approximation 
breaks down at low frequencies and large magnitude of M, [Fig. 3(b)]. 

Here we discuss some previous works using similar theoretical methods. Bleaney^° used 
the similar spin Hamiltonian to Eq. (1) in a strong field limit to examine the hyperfine struc- 
ture in paramagnetic salts. Friedman and Low calculated the resonant fields for Mn^+ doped 
into zinc fluosilicate (ZnSiF6:6H20) with the apphed field oriented at an angle with respect 
to the crystal axes, assuming that the Zeeman energy is much larger than the zero-field 
anisotropy energy. This approximation, which is the same as discussed in Appendix A of 
the present paper, is a proper assumption for Mn^+ in zinc fluosilicate since that system has 
a small crystal-field anisotropy. If we apply this assumption to the single-molecule magnet 
Mni2 however, we find that the resonant field within this approximation is in poor agree- 
ment with the exact diagonalization results. (See Fig. 4 and Appendix A.) The second-order 
corrections overestimate the exact results. For Mni2 at most intermediate frequencies, the 
Zeeman energy is not very large compared to the zero-field anisotropy, and the appropriate 
expressions for the resonant fields are our Eqs. (3) and (4). 

IV. DISTRIBUTION FUNCTIONS FOR THE RESONANT FIELDS 

Since we obtained the resonant fields as functions of g, D, and 9 in Eqs. (3) and (4), 
we now calculate numerically the probability-distribution functions (pdf) for the resonant 
fields, using Gaussian distributions in g, D, and 9, to predict the experimental EPR spectra. 
In our current study, the asymmetry in the spectra is concerned rather than the relative 
intensities of the spectra. The effect of the perturbation V due to the nonzero 9 (and the 
distribution in 9) on the relative intensities will be briefly discussed in Sec. V. The pdf of 
the resonant held, Fb{B), can be calculated as follows. 




(6) 
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fBMB,D,e)=fgM9*{B),D,9)/ 



(7) 



where B denotes the resonant field, Eq. (3) or (4), and g*{B) is obtained by solving the 
resonant field Eq. (3) or (4) for g?^ The function fg^D,e{g*{B), D,9) is the joint pdf of the 
three random variables g, D, and 9, calculated at g — g*{B). For simplicity, we assume 
that g, D, and 9 are statistically independent, so that the joint pdf factorizes. The double 
integration in Eq. (6) was performed numerically using Mathematica.^^ Figure 7 shows the 
pdfs of the resonant fields for a few transitions at low {v — 65 GHz and 66.135 GHz) and 
high {y = 200 GHz) frequencies. Here we use fixed values of the standard deviations of D 
and g, {ao — 0.02D and ag — O.OOSg, which are the same values used in Ref. 13) and two 
different values of the standard deviation of 9, uq fa 2.9° and 5.7°. 

As expected from the previous section (see Fig. 6), for = 2 ^ 1 [Figs. 7(a) and (b) 
for u = 66.135 GHz and u = 200 GHz, respectively], the long tail of each pdf appears on 
the right hand side of the maximum of the pdf (the peak field). For both frequencies, the 
long tails are recognizable even for ae ^ 2.9°, and they become prominent for ae fa 5.7°. 
For Ms — 3 ^ 2, at 1/ = 66.135 GHz, it is hard to visually recognize any asymmetry [Fig. 
7(c)] because the asymmetry is much smaller and the symmetric linewidth is larger than 
for Ms = 2 ^ 1. However, we see a small shift of the peak field towards higher fields as ae 
increases from 2.9° to 5.7°. As the frequency increases, the asymmetry effect is significantly 
enhanced [Fig. 7(d)]. For = — 4 — > —3, a possible asymmetry is expected at low 
frequencies only, as seen from Fig. 6(b). The peak field shifts slightly towards lower fields as 
a0 increases from 2.9° to 5.7° [Fig. 7(e)]. Other than that, it is hard to see any asymmetry 
in the pdf for the same reasons as for = 3 — > 2. For Mj — > + 1, as the magnitude of 
Ms increases, the asymmetry effect increases substantially [Compare Fig. 7(c) with 7(f)]. 

To quantify the asymmetry, we calculate the third central moment, {{B — (B))^), of 
the resonant-field distribution, and the difference, {{B) — -Bpeak), between the average field 
and the peak field, for the transitions shown in Fig. 7 (see Tables I and H). If a pdf is 
symmetric, its third central moment vanishes, and the peak field should coincide with the 
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average field. If a long tail appears on the right (left) hand side of the peak field, then the 
sign of the third central moment is positive (negative). For Mg — 2^1 and = 3 ^ 2, 
the signs of ((5 — (-B))^) and {{B) — -Bpeak) are positive, whereas for Ms — —4 — > —3 and 
Ms = — 10 ^ — 9 they are negative, in agreement with Fig. 7. 

V. DISCUSSION AND CONCLUSIONS 

As shown in Figs. 5-7, ag a; 5.7° is large enough to observe asymmetric tails in the 
EPR spectra for some frequencies and energy levels, but small enough to be realistic in 
experimental samples. We separately discuss the asymmetry in the spectra at a particular 
value of ag for Ms ^ Ms - 1 and for ^ + 1. For ^ - 1 (in the lower- 
energy potential well) , the asymmetry in the spectra is more pronounced for smaller Ms and 
higher frequencies, until the frequency becomes very high. At very high frequencies (about 
500 GHz or higher), the asymmetry in the resonant field increases weakly with increasing 
Ms, but the asymmetry in the spectra for large Mg would be masked by large symmetric 
linewidths. For Mg — > Ms + 1 (in the higher well), the asymmetry in the spectra is strong 
for large magnitude of Ms and low frequencies, because of high resonant fields. But for very 
low frequencies (below about 50 GHz) a different approximation scheme is needed than Eqs. 
(3) and (4). Since for high frequencies the EPR excitation energy becomes larger than the 
difference in adjacent energy levels, the transitions Ms — >■ Ms + 1 have narrow frequency 
or energy-level windows, in which asymmetric spectra are observable. Thus intermediate 
frequencies (between about 50 and 100 GHz) would be optimum for observation of the two 
different types of asymmetry in the spectra. 

So far we have concentrated on the asymmetry in the EPR line shapes without mention- 
ing intensities of the spectra. Here we briefly discuss the effect of the distribution in the 
tilt angle 9 on the intensities of the EPR spectra. At fixed values of D, g, and 6, the power 
absorbed between the Ms and the Ms — 1 energy level is written as 

dS 

= ^Ms\PMs,Ms] + ^Ms-i[Pm,-i,m,-i], (8) 
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where p{t) is the density matrix of the spin system, p{t)m'm — [pms,mJ is the 

change with time of the population in the M, level. To obtain [f>Ms,Ms]i density matrix 
equation^'' of the spin system [Eq. (1)] is used. For transition rates between different energy 
levels, the spin-phonon couphng Hamiltonian^^ is used. For more technical details, refer 
to Ref. 13, 34. Near resonance (angular frequency of the external transverse oscillating 
magnetic field UJ ~ {^Ms ~ ^Ms-i)/^), the power absorption becomes^^ 

dS_ ^ (^M,^LZ^|(M,|i?,5.|M, - 1)1^ A(i?,) {pM.Ms - PM.-i,Ms-i) , (9) 

where A(5^) is a Lorentzian line-shape function^^ and Pm^Ms population of the Ms 

level. The average power absorption can be calculated by averaging Eq. (9) over the Gaussian 
distributions m. and 6. The relative intensity of the average power absorption at a fixed 
resonance frequency is determined by the transition probability, the line-shape function, the 
population difference (Boltzmann factors), and spreads in D, and 6. The transition 
probability changes with the energy levels involved with the EPR transition, but it does 
not change due to the perturbation V caused by the nonzero 0. The reason is that as 
far as we are interested in the phenomena near resonance, only perturbation terms which 
oscillate with the frequency close to the resonance frequency are important to the transition 
probability. Thus, the interaction, Q{t) = BxSx{e^'^^ + e~'''^^)/2, between the spin system and 
the oscillating transverse magnetic field does contribute to the transition probabilities, 
but not the perturbation terms V. Then, the maximum height of the line-shape function 
A{Bz) is determined by a linewidth due to the finite lifetime of an excited state. The lifetime 
of an excited state does not change much with the perturbation terms V, so that the height 
is not much affected by the nonzero 9. The populations of the Mg and Mg — 1 level are 
not affected by the nonzero 6, because the populations are determined by the Boltzmann 
factor and small changes in the energy due to V are not noticeable in the Boltzmann factor. 
Therefore, we conclude that the relative intensity for an EPR transiton between specific 
energy levels is not affected by the perturbation V, and that the effect of the distribution 
in 9 on the relative intensity is solely due to the spread in 9. 
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To observe asymmetric line shapes in EPR experiments, one needs to optimize the exper- 
imental conditions. First, the asymmetry is more prominent for smaller Mg for Mg — > Mg — 1 
(or larger magnitude of Mg for Mg — > Mg + 1), but for those transitions the spectral inten- 
sity is generally poor at low temperatures because of the small populations in the excited 
states. Thus the sample temperature must be optimized to achieve a reasonable intensity. 
Second, one must avoid the level crossings at which levels in the two potential energy wells 
coincide.^^ If the EPR transitions happen to occur near level-crossing points, the spectra 
could pick up a large extra line broadening which could prevent one from observing the 
small asymmetry. As the frequency increases, adjacent energy levels move farther apart, 
so it becomes easier to avoid level-crossing points in the EPR transitions. Third, one must 
avoid asymmetries caused by experimental artifacts, such as the presence of standing waves 
in the EPR probe. ^'^^ Standing waves can cause severe line-shape distortion due to mixing 
between a dispersion spectrum and an absorption spectrum. This can be avoided by using 
the resonance cavity EPR technique. ^^^^ Fourth, the choice of frequency is important. To 
observe two different types of asymmetry (one in the direction of increasing field and the 
other in the direction of decreasing field) at a single frequency, the frequency must be lower 
than about 100 GHz. 

If the asymmetries are observed in spectra with the optimum experimental conditions, 
one can estimate the distributions in D and g from the linewidths of almost symmetric 
spectra. Then from the spectra with significant asymmetries, we can estimate how broadly 
the easy axes of the molecules are distributed in the sample. The width of the distribution 
in 9 may provide information on the defect concentration. Additionally the distribution in 
9 may partially contribute to the distribution in the tunnel splittings, which was recently 
proposed^ and measured experimentally,^^ because nonzero 9 produces the transverse terms 
relevant to tunneling from one potential well to another. 

In summary, we have theoretically examined the effect of defects on the line shapes of 
EPR spectra for field sweeps with the quasi-static external field along the crystal c axis at 
fixed frequencies in the range of 50 GHz to 500 GHz for a single crystal of the single-molecule 

14 



magnet Mni2. Static, random rotations of the magnetization easy axes of the molecules, due 
to defects, lead to asymmetries in the spectra. The strength and direction of the asymmetry 
depend on the frequency, the energy level, and the width of the distribution in the tilt angle 
9 (which depends on the defect concentration). With carefully chosen frequencies, energy 
levels, and the appropriate experimental technique, the asymmetry effect should be observ- 
able in future experiments. Comparison of the observed and calculated spectra could yield 
direct and quantitative information on the distribution in orientations and the concentration 
of defects in the samples. This information may also be useful to understand the origin of 
magnetization tunneling in the molecular magnet Mni2. Experimental examination of our 
predictions should be quite fruitful. 
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APPENDIX A 



The spin Hamiltonian shown in this Appendix is equivalent to the spin Hamiltonian (1), 
but the approximation used in the calculations are different from those used in Sec. III. 

To rewrite this spin Hamiltonian in terms of the spin operators along the crystal axes, 
we use the following rotation matrix, 



Sx I cos ip cos — cos 9 sin sin cos ip sin + cos 9 cos sin ip sin ip sin 9 

S„ = — sin -0 cos — cos 9 sin cos ip — sin -0 sin + cos 9 cos cos ip cos ip sin 9 
sin ^ sin — sin 9 cos cos 9 

where {9, 0, -0) are the three Euler angles. Then we obtain 



V 
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n = -D[sin^ e sin^ 5f + sin^ 9 cos^ Sl + cos^ ^ _ gj^s ^ gin cos -S^} 

— sin cos 6* cos 4> {Sbi Sc} + sin cos 6* sin {Sa, Sc}] — qiibBSc , (Al) 

where {A, B} is the anticommutator. 

To calculate the resonant fields, we here assume that the Zeeman energy is much larger 
than the zero-field anisotropy energy, so that the eigenvalues M'^ of Sc are good quantum 
numbers. Taking the Zeeman energy, —gfisBSc, as an unperturbed spin Hamiltonian (taking 
the quantization axis as the direction of the external magnetic field) and the other terms as 
small perturbations, we obtain the resonant fields as follows: 

S,es(M^ ^ - 1) 



The above results were also shown in Ref. 25. As shown in Fig. 4, the resonant fields in 
this approximation do not agree well with the exact diagonalization results, in contrast to 
Eqs. (3) and (4). For most intermediate frequencies, the Zeeman energy is not very large 
compared to the zero-field anisotropy energy, and the approximation of Eqs. (3) and (4) is 
therefore to be preferred. 



1 2M' — 1 



gfj^B^ 2 

D^sin^^cos^^ 



(431 + 24M^ - 24Mf ) 



2hi' 

^^^(-217-6M; + 6Mf)], 
1 . . 2M' + 1 



-[-hu- 



gi^B^ 2 

sin^ 9 cos^ 9 

H — 



-0(3008^9-1) 



+ 



2hv 
D'^ sin^ 9 



(-431 + 24M^ + 24Mf ) 



Shu 



(217 - 6M^ - 6Mf )] . 



(A2) 



(A3) 
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TABLES 

TABLE I. Third central moments of the resonant-field distribution, and the difference between 

the average and peak resonant fields at u = 65 GHz, 66.135 GHz, and 200 GHz for ag 2.9°, 
an ~ 0.02D, and ag « O.OO85 (see the solid curves in Fig. 7). The third root is taken for the third 
central moments to give the same units for both measures of the asymmetry. For Mg = —4, the 
asymmetry is too small to obtain a reliable moment and difference. 



frequency (GHz) 


energy level Mg 


^{{B - {B)r) (tesla) 


(B) - Bpeak (tesla) 


66.135 


2 


0.023 


3.02 xlO-3 


66.135 


3 


0.0052 


1.90 xlO"^ 


65 


-10 


-0.27 


-6.68 X 10-2 


200 


2 


0.030 


3.07 xlO-^ 


200 


3 


0.025 


2.40 xlO-3 


TABLE n. Third central moments of the resonant-field distribution, and the difference between 
the average and peak resonant fields at = 65 GHz, 66.135 GHz, and 200 GHz with the same 
values of an and ag as in Table I, but with ag = 5.7° (see dashed curves in Fig. 7). 


frequency (GHz) 


energy level Mg 


^{{B-{B)r) (tesla) 


(B) - Speak (tesla) 


66.135 


2 


0.068 


9.89 xlO-3 


66.135 


3 


0.016 


1.11 xlO-^ 


65 


-10 


-0.93 


-4.46 X 10-1 


200 


2 


0.098 


2.69 xlO-2 


200 


3 


0.091 


2.13 xlO-2 
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FIGURES 




FIG. 1. Schematic diagram of the magnetic core of the Mni2 molecule. Each molecule includes 
a tetrahedron of four Mn'^~^{S = 3/2) ions at the center, surrounded by eight Mn^"^(S' = 2) ions, 
yielding an effective ground-state spin of 5 = 10 (5 = 8x2 — 4x 3/2). Each molecule has a 
tetragonal symmetry. 
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FIG. 2. Schematic diagram of the potential-energy wells for the molecular magnet Mni2, which 
has an effective spin S = 10. Here Mg is the eigenvalue of the spin operator S^- For the lower energy 
well (the right-hand well), the difference between adjacent energy levels is g/isB + D(2Ms — 1). For 
the higher energy well (the left-hand well), the difference is —g^B^ + D{2Ms — 1). As the applied 
field increases, the difference between the energy levels Mg = —10 and Mg = +10 increases. 




-0.10 I ' — ' ' ' ' ' ' ' ' 1 

1 2 3 4 5 

Frequency (100 GHz) 




0.0 0.5 1.0 1.5 2.0 



Frequency (100 GHz) 

FIG. 3. The ratio between the second-order and the zero-order terms in the resonant fields, 
y(z^, Ms, 6), vs frequency at ^ 5.7° for the transitions (a) Mg ^ Mg - 1 and (b) Mg ^ Mg + 1. 
See details in the text. Equations (3) and (4) are valid when \Y{u,Mg,9)\ <C 1. In (a) Mg = 1 
represents the transition Mg = 1^0. In (b) Mg = —3 represents the transition Mg = — 3 ^ —2. 
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•— • Exact 




<^^^2nd-order in D 






^^^2nd-order in sin6 











0.02 0.04 0.06 0.08 0.1 

G (radian) 



FIG. 4. Resonant fields vs tilt angle ^ at =66.135 GHz for the transitions (a) Mg = 3 — > 2 and 
Ms = — 4 ^ —3, and (b) Mg = 2^1. The filled circles denote the exact results from numerical 
diagonalization of the spin Hamiltonian. The triangles denote the results from Eqs. (3) and (4). 
The diamonds denote the results from Eqs. (A2) and (A3). Eqs. (3) and (4) show good agreement 
with the exact results, in contrast to Eqs. (A2) and (A3). 
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0.15 




Energy level Ms 

FIG. 5. Asymmetry A{v,Ms,9 = 5.7°) vs Mg at various frequencies for the transitions (a) 
Ms-1 and (b) Mg ^ Mg + 1. See details in the text. 
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0.5 1 1.5 2 2.5 

Frequency (100 GHz) 

FIG. 6. Asymmetry A{u, Ms,9 = 5.7°) vs the frequency u for the transitions (a) Mg 
and (b) Mg Mg + 1. See details in the text. 
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FIG. 7. The probability distribution functions for the resonant fields for the transitions: (a) 
Ms = 2 ^ I, {c) Ms = 3 ^ 2, (e) M, = -4 ^ -3 at = 66.135 GHz, (b) M, = 2 ^ 1, (d) 
Ms=3^ 2 at u = 200 GHz, and (f) Ms = -10 ^ -9 at u = 65 GHz, when B\\c. Here we use 
the standard deviations of D as 0.02D, the standard deviation of g as O.OOSg, and the standard 
devation of 6 as 2.9° (solid curves) and 5.7° (dashed curves). 
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